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Quantum interference in coherent transport through single molecular rings may provide a mechanism to 
control current in molecular electronics. We investigate its applicability by using a single-particle Green function 
method combined with ab initio electronic structure calculations. We find that the quantum interference effect 
(QIE) depends strongly on the interaction between molecular tt states and contact a states. It is masked by a 
tunneling in small molecular rings with Au leads, such as benzene, due to strong tt-ct hybridization, while it is 
preserved in large rings, such as [18]annulene, which then could be used to realize QIE transistors. 



The use of single molecules as functional devices is the 
ultimate end of the ongoing trend toward miniaturization of 
electronic circuits 1,2,3 . Despite significant progress made in 
the last decade, several issues still challenge the realization of 
molecular electronics, notably the sensitivity and control of 
molecule-lead contacts 4,5,6 . Recently Cardamone et al. 7,8 pro- 
posed a novel mechanism to control electron transport through 
single molecular rings: the current is determined by the degree 
of destructive or constructive quantum interference between 
the two paths around a symmetric molecule. Such interference 
can be controlled by a third terminal providing elastic scatter- 
ing or dephasing. Their use of a quantum interference effect 
(QIE) in the completely coherent quantum regime builds on 
previous work in semiconductor nanostructures 9, 10 and molec- 
ular nanostructures 11,12,13,14,15,16 . Since the QIE stems essen- 
tially from the symmetry of a molecular device, it should not 
be affected significantly by the structure of the molecule-lead 
contact. We present the first study of this novel mechanism in 
realistic ab initio calculations. 

QIE controlled molecular electronics was proposed 7 ^ us- 
ing a model calculation for a benzene ring in which only the 
7r molecular states were considered. The connection between 
the leads and the molecular states was treated phenomenolog- 
ically, and the possibility of induced structural relaxation was 
ignored. As these simplifications could have a big effect on 
the conductance, it is critical to investigate the validity of the 
QIE mechanism in realistic systems and to evaluate the factors 
affecting it. 

In this paper, we investigate the QIE mechanism by per- 
forming quantum transport calculations for two molecular 
rings, one small - benzene - and one large - [1 8]annulene. We 
adopt a standard single-particle Green function methoctfJl 
combined with ab initio electronic structure calculations in 
which all the above-mentioned factors are fully taken into 
account. Two different electrodes (leads) are studied: gold, 
which has ir-a molecule-lead coupling, and a metallic (5,5) 
carbon nanotube (CNT), which has strong ir-ir coupling. 

Our calculations show that the survival of the QIE mecha- 
nism depends strongly on the interaction between the molec- 
ular 7r states and the contact a states. It is masked by a tun- 
neling in the small benzene/ Au system because of the strong 
tt-ct hybridization, while it is preserved in the large [18]annu- 
lene/Au system, which then could be used to realize QIE tran- 
sistors. With the CNT leads, we find that strong ir-ir molecule- 
lead coupling can modify significantly the QIE but does not 



destroy it. 

Our method for calculating electron transport through 
molecular junctions has been described in detail previously 18 ; 
it combines a Landauer approach to transport with ab initio 
electronic structure method s 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 . i n particular, 
we use here density functional theory (DFT) 27 in the local 
density approximation (LDA), hybrid DFT in the version of 
B3LYB2&2?-, and Hartree-Fock (HF) theory. Several energy 
functionals are used to make sure that the conclusions are gen- 
erally valid, because to date no functional is completely ac- 
curate for transport calculations. Usually DFT-LDA and HF 
are the two extremes, the former overestimates and the latter 
underestimates the conductance 30 . The wave functions are ex- 
panded using a Gaussian 6-31 1G** basis set for the C, H, and 
S atoms, and a CRENBS basis set for Au atoms 31 . 

Practically, the lead-molecule-lead system is divided into 
three parts: left lead, right lead, and device region. The lat- 
ter contains the molecule plus parts of the leads to accom- 
modate the molecule-lead interaction (see structures below). 
The self-consistent DFT or HF Hamiltonian of the device 
region plus the self-energies of the two semi-infinite leads, 
^l,r(E), are used to construct a single-particle Green func- 
tion, Gd(£), from which the transmission coefficient as a 
function of energy is calculated: T(E) = Ti[T L G D T R G^ D ] 
where T LyR (E) = i[H LtR (E) - ^ LjR (E)] is the coupling to 
the left or right lead. The conductance, G, then follows from 
a Landauer- type relation. 

For the benzene ring, we consider only gold leads. The self- 
energy due to the Au within the leads is treated in the wide 
band limit (WBL) approximation: T(E) = — 7I with 7 = 3.0 
eV— . But note that the key molecule-lead coupling is included 
explicitly in the Hamiltonian because 9 Au atoms are included 
in the central device region. Since the density of states of gold 
is quite flat around the Fermi energy, the use of the WBL for 
Au atoms is reasonable. For the large [18]annulene ring, we 
consider gold leads as for benzene as well as (5,5) metallic 
carbon nanotube leads for which the self-energy is obtained 
by ab initio calculation of atomic leads. For all systems the 
atomic structure of the junction, including the molecule-lead 
separation, is optimized 33 by minimizing the atomic forces on 
atoms to be smaller than 0.02 eV/A. 

Let us start with the small system - the benzene ring in 
Fig. 1(a). All of the Au atoms included in the device region 
are shown; the WBL self-energy is applied to these atoms. 



2 





FIG. 1: (color online) Optimized structures of the Au-S -benzene- S- 
Au system (device region only) with (a) the (1,4) constructive and 
(b) the (1,3) destructive interference configurations, (c) The alkane 
configuration used for comparison. (d),(e) Transmission functions 
obtained with different energy functionals. For the HF calculation 
we consider cases both with and without the lead-induced molecu- 
lar relaxation, as indicated. Note the absence of a clear difference 
between the (1,3) and (1,4) cases. 



The S atom is situated on a hollow site. Two structural config- 
urations are considered: (1,4) and (1,3) in Figs. 1(a) and 1(b), 
respectively. Configuration (1,4) has constructive interference 
in the all-7r model (phase difference between the two paths 
around the ring is 0), while (1,3) has destructive interference 
(phase difference is 7r) 7,8 . Still within the all-7r model, the 
non-equilibrium many-body physics was subsequently stud- 
ied in detail for these two configurations 34 . 

The transmission T(E) found using several functionals is 
shown in Figs. 1(d) and (e). First, note that the transmission 
gap depends on the functional as expected: it increases in 
the order LDA, B3LYP, and HF. Also, the LDA equilibrium 



FIG. 2: Local density of states (LDOS) for the (1,3) benzene/ Au 
system using LDA. (a) Top view and (b) side view for the energy 
window [-2.5, -1.0] eV (HOMO resonance) in Fig. 1 (d). (c) Top 
view for the energy window [-0.2, 0.2] eV. The influence of the a 
bonds is clear in (a) and (b). 



conductance for the (1,4) configuration is consistent with pre- 
vious calculations 35 ! 36 ! 37 , about 0.1 Go. The key result here 
is the comparison between the (1,4) and the (1,3) configura- 
tions. The all-7r model predicts that the (1,3) configuration 
has a transmission node at the Fermi energy due to total de- 
structive interference 7,8 . However, in our more realistic cal- 
culation, this transmission node does not occur. For LDA and 
B3LYP G(i,3) is smaller than G(i 4 ) by a factor of about 2; 
for HF the order is even reversed, G(i,3) > G{i,a) • 

Several factors may destroy the perfect destructive interfer- 
ence present in the model of Refs. and|H: (i) lead-induced 
structural relaxation which breaks the symmetry, (ii) the influ- 
ence of the a states, and (iii) beyond nearest-neighbor interac- 
tions which cause other paths. Let us first look at the effect of 
structural relaxation using the HF calculation: results for re- 
laxed and unrelaxed structures are compared in Fig. 1(e). The 
effect on T(E) is clearly minor. Because of the symmetry of 
the additional paths, the effect of (iii) should also be small. 
Thus, we conclude that factors (i) and (iii) are not important. 

The first hint of the strong effect of a hybridization comes 
from the difference between the transport gap in Fig. 1 and 
the HOMO-LUMO gap of an isolated benzene molecule. In 
LDA, for instance, the former is ~4.1 eV while the latter is 
5.2 eV. This suggests that the HOMO and LUMO resonances 
may not be dominated by benzene tt states. To check this, we 
plot in Fig. 2 the LDA local density of states (LDOS) both for 
the HOMO resonance (— 2.5 eV to —1.0 eV) and for near the 
Fermi energy (—0.2 eV to 0.2 eV). Figs. 2(a) and (b) show 
that the HOMO resonance has large contributions from the 
contact Au-S states. Even in the benzene ring the a character 
is remarkable. Because of this strong ir-a hybridization, trans- 
port near Ep consists largely of tunneling through the a states 
rather than through the 7r-ring, as clearly shown by the LDOS 
around Ep in Fig. 2(c). To further examine this conclusion, 
we replace the benzene ring with a a-bonded C3H6 molecule 
[Fig. 1(c)]. T(E) in Fig. 1(d) (green line) shows that the con- 
ductance is very close to that of the (1,3) system, suggesting 
a similar transport mechanism. This shows directly that the 
conductance in the (1,3) configuration mainly comes from a 
tunneling which masks the QIE effect. 

Next we turn to the larger [18]annulene ring shown in Fig. 
3. We consider two configurations which yield destructive 
interference in the 7r-only model 7,8 , (1,7) and (1,9), and one 
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FIG. 3: (color online) Optimized structures for Au-S-[18]annulene- 
S-Au (device region only) with (a) (1,7), (b) (1,9), and (c) (1,10) 
lead configurations, (d) The alkane structure used for comparison. 
Transmission functions are given in (e) LDA, (f) B3LYP, and (g) HF. 
Note the clear difference between configurations with destructive and 
constructive interference. 



constructive configuration, (1,10). From the T(E) given by 
LDA, B3LYP, and HF one sees the same trend in the transport 
gap and in the equilibrium conductance: the smallest conduc- 
tance is given by HF with the largest gap, which is about one 
order of magnitude smaller than the LDA result, in agreement 
with a previous result for other long-chain molecules 30 . The 
key feature here is that the conductance for (1,7) and (1,9) is 




(a) 
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FIG. 4: LDOS for the (1,7) [18]annulene/Au system from LDA. 
(a) Top view and (b) side view for the energy window [-1.0, -0.2] eV 
(HOMO resonance) in Fig. 3 (e). (c) Top view for the energy window 
[-0.2, 0.2] eV (Fermi energy). The 7r-orbital signature is clear. 



much smaller than that for the (1,10), by about two orders of 
magnitude. This indicates that, in contrast to benzene, the 
ideal destructive interference is largely preserved. 

To understand this, we show in Fig. 4 the LDA LDOS near 
both the HOMO resonance (-1.0 eV to -0.2 eV) and the Fermi 
energy (-0.2 eV to 0.2 eV). One sees that the HOMO reso- 
nance is dominated by the [18]annulene tt states, and trans- 
port near the Fermi energy is basically through the 7r-ring. 
Because of the near perfect destructive interference of the tt 
transmission 7,8 , the residual conductance should come from 
a-bond tunneling. To check this, we replace the [18]annulene 
ring with a a-bonded C7H14 chain [see Fig. 3(d)]. The T(E) 
[green line in Fig. 3(e)] yields a conductance very close to that 
of (1,7), confirming our view. 

How to understand the difference in behavior between the 
small benzene ring and the large [18]annulene ring? The an- 
swer lies in the energies of the molecular levels. In LDA, for 
example, the energy of benzene tt states (HOMO) is -6.5 eV 
which is in the energy window of the Au-S bonding at the 
contact, from -7.0 to -5.7 eV. The larger size of [18]annulene 
means that its tt states (HOMO) are higher in energy, at -5.0 
eV, which is well above the energy window of the Au-S bond- 
ing. Consequently, in [18]annulene, the tt orbital becomes the 
frontier orbital which therefore dominates transport. 

From this picture one can see that the advantage of using 
larger conjugated molecules is two-fold: First, the resulting 
7r-dominated transport preserves the simple quantum interfer- 
ence effects. Second, the longer length suppresses a tunnel- 
ing, leading to a larger on/off ratio. 

While Au electrodes have some advantages, they have a se- 
vere disadvantage as well: the contact transparency is poor 
(despite the strong chemical bond), leading to a small equi- 
librium conductance even in the case of constructive interfer- 
ence. For a more favorable case with strong coupling near the 
Fermi energy, we use a metallic (5,5) carbon nanotube as a 
lead, connecting to the molecule through a 5 -member-ring 38 . 
The optimized structures are shown in Fig. 5 for the (1,10) and 
(1,7) configurations; both have a co-planar structure which 
provides overall conjugation and thus strong molecule-lead tt- 
ix coupling. 

The strong molecule-lead coupling leads to a metal-induced 
gap state traveling into the molecule 38 . For short conjugated 
molecules, the two states originating from the two leads will 
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FIG. 5: Optimized structures of the 5 -member-ring-connected 
(5,5)CNT-[18]annulene-(5,5)CNT junctions with the lowest total 
energies: (a) constructive (1,10) and (b) destructive (1,7) configura- 
tions. Their transmission functions for different molecule-lead sepa- 
rations (indicated in the legend) are given in (c) and (d), respectively. 
Note the resonance peak around the Fermi energy in (c) and the anti- 
resonance peak in (d). 



meet in the molecule, causing a resonance peak in T{E) 
around E F . This is just what we see in Fig. 5(c) for the 
constructive (1,10) configuration. [For T(E) here, we do 
not use the WBL, but rather find the lead self-energy from 
an atomistic ab initio calculation.] The resonance peak de- 
cays as we increase artificially the molecule-lead separation 
(weaken the coupling), as shown in Fig. 5(c). For the de- 
structive (1,7) configuration [Fig. 5(d)], an anti-resonance 
peak appears at E F , indicating that the destructive interfer- 
ence applies in the strong coupling limit to the metal-induced 
gap state (though the anti-resonance is rather narrow). For 
larger molecule-lead separation (weaker coupling), the anti- 
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FIG. 6: (color online) (a) A schematic drawing of a [18]annulene- 
based field-effect QIE transistor. The source-drain current is con- 
trolled by the gate voltage V g which causes elastic scattering in the 
gated region. The transmission functions for zero bias (Vb=0) and 
different V g are given in (b). 

resonance becomes wider and deeper; thus, the molecule- 
lead coupling affects significantly the quantum interference 
but does not destroy it. Note that here the metal-induced gap 
state leads to a much larger constructive conductance and a 
larger on/off ratio compared with the case of Au leads. 

The effective destructive quantum interference can be used 
to control electric current through single molecular devices. 
To demonstrate this explicitly, we propose a model device 
structure, shown in Fig. 6, using a Au-S-[18]annulene-S-Au 
(1,7) junction, where a local gate potential (due to an STM tip, 
for instance) is applied to part of the molecular ring, which 
will destroy partially the destructive interference because of 
the additional scattering, depending on the gate voltage, V g . 
The calculation is carried out using the B3LYP functional. 
The computational techniques for transport in the presence of 
a local gate potential shift has been described previously 39 . 
Fig. 6(b) shows the transmission functions for different gate 
voltages. One can see that a small gate voltage can mod- 
ify significantly the transmission function of the device, both 
shifting and greatly distorting the transmission valley. Most 
importantly, a small gate voltage can increase the equilibrium 
conductance by up to three orders of magnitude - a good be- 
havior for field-effect QIE transistors. 

We would like to finish this Letter with a comment 
about the effect of molecular thermal vibration on QIE. 
For non-interference-related electron transport through sin- 
gle molecules which are not very long, molecular vibration 
should have a very small effect on the I-V characteristics be- 
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cause the transport is basically ballistic as long as extemely 
soft vibrational modes (which may lead to comformational 
changes) are absent. For QIE, the relevant vibrational modes 
are those which change the length difference between the two 
paths. These modes therefore involve changing the C-C bond 
length and have very high frequencies. As a result, the effect 



of molecular vibration on QIE for room temperature can be 
expected to be small. 

We thank Charles Stafford for a series of valuable conver- 
sations. This work was supported in part by the NSF (DMR- 
0506953). 
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